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Abstract. Following Manin's approach to renormalization in the theory of computation, we 
investigate Dyson— Schwinger equations on Hopf algebras, operads and properads of flow charts, 
as a way of encoding self-similarity structures in the theory of algorithms computing primitive 
and partial recursive functions and in the halting problem. 



1. Introduction 

There are many deep connections between theoretical physics and information theory, and in par- 
ticular the theory of computation, see for instance the account given in [T] . In the recent papers pT| 
and [22] , Manin developed a new approach to the theory of computation and the halting problem, 
based on importing ideas and techniques from the Hopf-algebraic formulation of renormalization 
in perturbative quantum field theory, |3], [9], [H]. The purpose of this paper is to show that the 
Hopf algebra of flow charts introduced by Manin in [21] exhibits self-similarity structures given by 
solutions of combinatorial Dyson-Schwinger equations, defined as in perturbative renormalization, 
[2] ; [10] ; [H] : [H] , [IS] ; [30] • These Can be thought of as a notion of "equations of motion" in the 
theory of computation. 

Dyson-Schwinger equations [7], |26j are a formulation of equations of motion in perturbative 
quantum field theory, expressed in the form of relations between Green functions. In recent years, 
an algebraic formulation of the combinatorial structure of perturbative renormalization for scalar 
field theories was developed, starting with the work of Kreimer on the Hopf-algebraic structure of 
renormalization TF , followed by the Connes-Kreimer formulation p| of the BPHZ renormalization 
procedure, and the formulation in terms of a Riemann-Hilbert correspondence for categories of 
differential systems in Connes-Marcolli [5] , [6] . The Connes-Kreimer formulation of renormalization 
was given a very general algebraic formulation in terms of Rota-Baxter algebras in the work of 
Ebrahimi-Fard, Guo and Kreimer [9]. Correspondingly, the Dyson-Schwinger equations were also 
given a combinatorial form, reflecting the Hopf-algebraic structure of renormalization, in the work 
of Bergbauer, Kreimer, and Yeats, [2], [19], [30] and more recently with an extensive study by 

Foissy, [n], m- 

In this paper we investigate the formulation of Dyson-Schwinger equations in the context of 
Manin's approach to renormalization and computation. In ^2] we review some variants of the 
construction of a Hopf algebra of flow charts. These are diagrams, consisting of decorated planar 
rooted trees, that compute primitive recursive functions. In ij3] we discuss grafting operators and 
Dyson-Schwinger equations for flow charts, as a way of identifying certain self-similarity structures 
in the computation of primitive recursive functions. We review some known results about existence 
and uniqueness of solutions and on conditions under which the coordinates of the solution generate 
a Hopf subalgebra. We show that one can also consider ideals generated by the coordinates of 
solutions. Under the condition that the grafting operator is a cocycle, these can give rise to Hopf 
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ideals. The quotient Hopf algebra determines a subgroup of the affine group scheme dual to the 
(commutative) Hopf algebra of flow charts. This subgroup defines a Galois group in the sense 
of Yanofsky's Galois theory of algorithm [29] , consisting of those symmetries that are compatible 
with the self-similarity structure imposed by the Dyson-Schwinger equation. We then show an 
explicit example of a Dyson-Schwinger equation based on computation by binary trees that does 
not satisfy the Hopf algebra condition. In order to accommodate this type of example, that is 
relevant to the setting of computation, we show that it is convenient to take an operadic viewpoint 
on Dyson-Schwinger equations, which was already suggested in the work of Bergbauger-Kreimer 
[5]. In 21 we show that Dyson-Schwinger equations can be formulated in the operad setting, using 
a family of grafting operations, which is not required to satisfy a cocycle condition. We consider 
Dyson-Schwinger equations on an operad of flow charts. We prove the existence and uniqueness 
of solutions for these operadic Dyson-Schwinger equations and also for a properad version. The 
extension from operads to properads is motivated by the more recent formulation, given by Manin 
in [24], of a properad of flow charts based on directed (acyclic) graphs instead of rooted trees. 
Finally, in !j5] we discuss the algebraic Feynman rule associated to the halting problem, as in [22] 
and its BPHZ renormalization and the role of Dyson-Schwinger equations in the halting problem. 

2. Primitive recursive functions and the Hope algebra of flow charts 

We recall here the setup of Manin i21j on the Hopf algebra of flow charts. First recall from [20] 
the following facts about primitive recursive functions. 

2.1. Primitive recursive functions. As in §V.2 of [20], we consider the class oi primitive recursive 
functions as generated by the basic functions 

• Successor s : N — > N, s{x) ~ x + 1; 

• Constant c" : N" ^ N, c"{x) = 1 (for n > 0); 

• Projection 7rf : N" ^- N, tt'/{x) ^ x, (for n > 1); 

With the elementary operations 

• Composition (substitution) C(„,^,„,p): for / : N" ^ N", g : N" ^ W, 

gof:W-^W, V{g o f) ^ f-\V{g)y, 

• Bracketing (juxtaposition) b(fe. for fi : N™ N"% i = 1, . . . , fc, 

/ - (A, . . . , /fc) : N" ^ N"l+■■■+"^ V{f) = P(/i) n • ■ • n Vifk); 

• Recursion c„: for / : N" ^ N and g : N"+2 N, 

/z(xi, . . . ,3^71, 1) . f {x\^ . . . ,X^), 

h{xi, . . . , x„, fc + 1) := g{xi, .. .,x„,k, h{xi, . . . , x„, fc)), fc > 1, 

where recursively (xi, . . . , a;„, 1) G T>{h) iff (a:i, . . . , x„) G T^if) and {xi, . . . , Xn, fc + 1) G 
2?(/i) iff {xi,...,Xn,k,h{xi,...,Xn,k) G X'(g). 

2.2. Hopf algebras of decorated rooted trees. We recall here the construction of Hopf algebras 
of decorated rooted trees and of planar decorated rooted trees and their main properties, see [11] 
and references therein. Throughout this paper vector spaces and algebras are over a field IK of 
characteristic zero. 
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2.2.1. Rooted trees. A rooted tree is a finite graph r, whose geometric reahzation |t| is simply 
connected, defined by combinatorial data (Fr, Vr,Vr, Srjjr), with Fr a set of half-edges (also called 
flags or tails) , Vr the set of vertices with a distinguished element Vr € Vt (the root) , a boundary 
map dr : Ft that associates to half-edge its boundary vertex and an involution jV - Ft — > -Ft , 

— 1 that performs the matching of half-edges that determines the edges of t. The resulting 
graph has edges (or internal edges) given by pairs of half-edges matched by the involution, and tails 
(or external edges) given by half-edges that are fixed by the involution. We denote by Et the set 
of internal edges and by E^^* the set of external edges. 

2.2.2. Orientations. We consider rooted trees as endowed with the orientation that the root vertex 
as the output, namely where all edges are oriented in the direction of the unique path to the root. 

2.2.3. External edges. We also assume that the root vertex has an outgoing half-edge and a number 
of incoming edges, while all other vertices have one outgoing edge and a number of incoming edges 
and possibly a number of incoming tails. Each leaf of the tree has an incoming tail. 

2.2.4. Planarity. A planar rooted tree is a rooted tree r together with a fixed embedding t : |t| 
of its geometric realization in the plane. 

2.2.5. Decorations. A (planar) decorated rooted tree is rooted tree and with a map (pv '■ Vr Vy 
of the set of vertices to a set Vy of vertex-labels and a map <j)F ■ Fr Vp to a set Vp of labels of 
flags. The assignment of a decoration {(j)v, (j^p) to a rooted tree r is in general subject to constraints. 
For example, if / and /' in Fr are matched by the involution jr to form an edge e £ Er, then they 
should carry the same decoration (j)F{f) = (I^F{f')- This defines a labeling of edges. We will see 
other constraints below, in the specific case of the Hopf algebra of flow charts. We also consider 
the forgetful map that forgets the flag labels Vp and retains the vertex labels Vy and a further 
forgetful map to unlabeled trees. 

2.2.6. Noncommutative Hopf algebra of planar rooted trees. Given a (planar) rooted tree r, an 
admissible cut C of r is a modification of the involution jr that cuts a subset of internal edges 
Bi G Er into two flags fi, f[. Namely, instead of having jr{fi) = fi, one modifles the involution so 
that the new j(/i) = fi and i{fi) = f', where the set of edges is chosen in such a way that every 
oriented path in r from a leaf vertex to the root contains at most one Cj. The new graph obtained 
using the involution j is a forest 



consisting of one (planar) rooted tree pc{'^) containing the root vertex of t and a finite disjoint 
union of other oriented planar trees TTc{r) — Ui7rc,i(T), where each tree TTc,i{T) has a single output 
external edge, to which we assign as root the boundary vertex of this output. 

The noncommutative Hopf algebra of planar rooted trees "H""^ is defined, as an algebra, as the free 
algebra generated by the planar rooted trees, with the coalgebra structure given by the admissible- 
cuts coproduct 



c 

The coproduct is coassociative and the Hopf algebra is graded by the number of vertices of rooted 
trees, hence the antipode is defined inductively by the formula 



C{T)=pc{T)Unc{r) 



(2.1) 
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with a;', x" terms of lower degree. See [13], [H] for more details on this Hopf algebra. An element 
a; in a Hopf algebra is primitive if A (a;) — l®x + x®\. 

2.2.1. Commutative Hopf algebra of rooted trees. In the Hopf algebraic approach to perturbative 
renormalization, as formulated by Kreimer |15j and Connes-Kreimer f3| , one considers a commuta- 
tive Hopf algebra of rooted trees decorated by Feynman graphs, with the admissible-cuts coproduct 
(|2.ip . or else a commutative Hopf algebra of Feynman graphs, with a coproduct that corresponds 
to inclusions of subgraphs (subdivergences) and quotient graphs. 

In the case of planar rooted trees, one can also consider a commutative Hopf algebra %, which 
is a quotient of H"'^ by imposing the commutativity of multiplication. The generators of H. are still 
planar rooted trees, but monomials in these generators no longer identify with embedded forests. 
The grading and antipode are as above, induced by those of H"'^. 

2.3. Manin's Hopf algebra of flow charts. We consider, as in [21], the set of planar labelled 
rooted trees, where the label set Vy of vertices is given by the set of elementary operations 

(composition, bracketing, recursion) and the label set "Dp of the flags is the set of primitive recursive 
functions. Notice that, because the inputs of the composition, bracketing and recursion operations 
are ordered, we need to work with planar rooted trees. 

The data [(f)v{v), 4>F{fi)) of a labelings (j)v[v) £ Vy and (t)F{fi) of a vertex v £ Vt and of all the 
half-edges /; G with d{fi) — v are admissible if they satisfy the following conditions: 

• If (jjviv) = C(m,n,p), then V must have valence three; the labels hi = (t>Fifi) and /i2 = 4'F{f2) 
of the two incoming flags must have domains and ranges hi : N™ N" and ft.2 : N" W 
and the outgoing flag must be labeled by the composition h2 hi ~ C(„i „ p)(/ii, ^12). 

• If (jiviv) = t„, then v must have valence three; the labels hi — 0f(/i) and ft.2 = 0f(/2) 
of the two incoming flags must have domains and ranges hi : N" N and /i2 : — N 
and the outgoing flag must be labeled by the recursion h — t„(/ii, h2). 

• If 4>v{v) = b(fc_m,ni) then v must have valence fc + 1 and all the functions hi — (f>Fifi) 
associated to the incoming flags must have domain in the same N™. The outgoing flag 
must be labeled by the bracketing / = (/i, . . . , Z^) = b(fc^„^„.) (/i, . . . , fk). 

Definition 2.1. An admissible labeling of a planar rooted tree t by primitive recursive functions 
is a collection of labelings 

{{(t)v{v),4'F{fi)) \v £Vr and £ Fr with d{f.,) = v} 

such that the admissibility condition above is satisfied at all vertices and such that 4>F{f) — <f>F{f') 
whenever jr{f) ~ f . 

A vertex labeling of a planar rooted tree r by elementary operations is a collection of labeling 
{0y (u) G {c, r, bfe} I V £ V-r} with the only constraint that labels of type c and r can only be assigned 
to vertices with two incoming flags and one outgoing flag and labels of type bk can be assigned to 
vertices with k incoming flags and one outgoing flag. 

Via the forgetful map that keeps the vertex labels and forgets the flag labels, a planar rooted tree 
with an admissible labeling by primitive recursive functions determines one with vertex labeling 
by elementary operations, while not all vertex labeling by elementary operations on a given planar 
rooted tree will admit a compatible flag labeling. Notice that, in the case of vertex labelings, we no 
longer distinguish between labels C(„i^„.p), with different choices of m, n and p, since we do not flx 
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the flag labels: all these vertex labels in Vy correspond to the same label c. Similarly, all the labels 
r„ of admissible labeling correspond to the same label r of vertex labeling, and all the b(^k,m,ni) 
correspond to the same bk- We retain the different bk as vertex labelings as those correspond to 
different valences of the vertex. 

Thus, a planar rooted tree with an admissible labeling by primitive recursive functions can be 
interpreted as a flow chart that computes the output function starting from the functions associated 
to the incoming external edges according to the operations performed at the vertices. 

A planar rooted tree r with a vertex labeling by elementary operations, on the other hand, 
should be thought of as a template of a possible computational scheme: when applied to a tuple 
of primitive recursive functions {hi, . . . , hk), with k the number of incoming flags of r, it either 
computes an output functions h, if the inputs (hi, . . . , hk) together with the assigned vertex labels 
determine an admissible labeling of r by partial recursive function, or else we set its output to be 
the empty function. 

We then define two slightly different versions of the Hopf algebra of flow charts, depending on 
the type of labeling that we want to use on trees. 

Definition 2.2. The noncommutative Hopf algebra 'H^^^ -p of flow charts is the free algebra gen- 
erated by the planar rooted trees with admissible labelings by primitive recursive functions, with 
the coproduct (j2.1[) . We also denote 6?/ Hfiow,'P the commutative quotient ofW^'^^^. Similarly, the 
noncommutative Hopf algebra "Hg^^ y is the free algebra generated by planar rooted trees with vertex 
labeling by elementary operations, again with the admissible- cuts coproduct. We denote by Hflow.v 
its commutative quotient. 

Notice that if r has an admissible labeling, then the rooted trees pc{t) and TTc,i{T) also inherit 
induced admissible labelings, so the Hopf algebras are well defined. There is a Hopf algebra homo- 
morphism 'H'^^^ p — >■ Hflow,V: determined by forgetting the fiag labels, but as observed above not 
all generators of Hflow.v are in the image of this map. Since the rooted trees in Hfiow,v have no 
labels attached to flags, we can represent them as rooted trees without external edges. This will be 
implicitly done in the description of the grafting operators in Lemma 13.11 below. 

In the following, whenever a statement applies to both types of tree decorations on the Hopf 
algebras of flow charts, we will drop the V and V subscripts. 

2.3.1. Binary versus k-ary operations. It is possible to reformulate the Hopf algebra of flow charts 
by using only binary trees. In fact, it is shown in §2 of [28 that the description of primitive recursive 
functions in terms of basic functions (successor, constants and projections) and elementary opera- 
tions (composition, bracketing and recursion), as recalled above, can be reformulated as repeated 
applications of binary operations of these same forms. In fact, the composition and recursion op- 
erations already are binary (they label vertices with two input fiags and one output) so one only 
needs to verify that the bracketing operation b(k^m,ni) assigning to k functions fi : N™ — > N"S 
i = 1, . . . , k their juxtaposition / = (/i, . . . , fk) : N™ — s> N"^^ is a composition 

b{k,m,ni) = &(2,m,ni,n2H hrifc) ° ' ' ' ° b(2,m,nfc_ i ,nfc ) ■ 

This replaces a single vertex with k inputs and one output with k vertices with two inputs and one 
output. Thus, without loss of generality in the algorithmic representation of primitive recursive 
functions, one can consider a sub- Hopf algebra Hg^^ C Hg^^ generated by binary planar rooted tree 
with an admissible labeling by primitive recursive functions, and the corresponding commutative 
quotient nl^^ C "Hs OW ■ 
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One the other hand, one can also extend the two binary operations of composition and recursion 
to fc-ary operations for arbitrary k > 2. In the case of composition, we can introduce a new set 
of vertex labels C(^k,ni) that correspond to the fc-ary compositions C(k,m,ni)(hi) = hk ° ■ ■ ■ o hi of 
functions hi : N"'"^ — >■ N"', for i — 1, . . . , k, with no = m. In the case of recursion, for fc > 1, one 
can consider (fc + l)-ary recursions where the recursion depends on k initial conditions. Namely, 
we define rfc+i^„ as the operation that computes the function h — Xk,n{hi, ■ ■ ■ , hk+i) : — s- N, 
from the input /i^ : N" -)> N for i = 1, . . . , fc and hk+i : -> N, as 

h{xi, . . . ,a;„, 1) = hi{xi, . . . ,a;„), 

(2-2) h{xi, ...,Xn,k) = . . . , x„), 

h{xi, . . . ,Xn,k + e) = hk+l{xi, ...,Xn,hi{xi,...,Xn),---,hk{xi,...,Xn),k + e-l), 

for £ > 1 

These fc-ary operations are clearly reducible to compositions of binary operations, hence any 
primitive recursive function that can be computed using these operations can also be computed by 
the presentation in terms of binary compositions and recursions. However, as we shall discuss more 
in detail below, it is convenient to include these explicitly in the structure of the Hopf algebra of 
flow charts in order to have well deflned grafting operators that give rise to Hochschild cocycles 
and to Dyson-Schwinger equations. 

Thus, we propose another minor modiflcation of Manin's Hopf algebra of flow charts, as follows. 

Definition 2.3. We denote by Hfl^^ -pi and "H^^^ y/ the noncommutative Hopf algebras of flow 
charts where the vertex label sets is enlarged to contain all the k-ary compositions and recursions 
C(fc,ni) and Vk,n, and the admissibility condition is restated accordingly. We denote by Hfiow.v and 
'Hflow,v' their commutative quotients. 



3. Dyson-Schwinger equations in the Hopf algebra of flow charts 

Recall here the formulation of combinatorial Dyson-Schwinger equations in Hopf algebras of 
decorated rooted trees, following [5], [33 and especially [TU], [H], [H]. We then focus on the 
speciflc case of the Hopf algebras of flow charts described above. 

3.1. Insertion operators for decorated planar rooted trees. We work here with the Hopf 
algebra of flow charts 'HJiq^ y/ . By analogy with the setting of perturbative renormalization [2] , 
[30] . we deflne, for each type of elementary operation 5 G {b, c, r}, a grafting operator B'^ deflned in 
the following way. Given a monomial T in Hgo^ y/ , which consists of a forest with vertex labeling 
by elementary operations, we deflne Bg'{T) as the sum of planar graphs obtained by adding a new 
root vertex with a number of incoming flags equal to the number of trees in T and a single output 
flag, with the vertex decorated by S. 

3.1.1. Hochschild 1-cocyles. In the case of perturbative renormalization, it is well known |? that 
the analogous grafting operators Bf , with vertex decorations 5, deflne Hochschild 1-cocycles. 
The 1-cocycle condition for the operator Bf then consists of the property that 

(3.1) /^Bj = {id®Bj)/^ + Bj ®l. 

Lemma 3.1. The operators B^ on TLg^^ y,, with S G {b, c,r}, satisfy the 1-cocycle condition p.ip . 
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Proof. This case works exactly as in the renormalization setting, see [2]: we reproduce the argument 
here for the reader's convenience. For an element x in a Hopf algebra H with coproduct 

A(a;) = l(K)a; + a;®l + ^a;'(8) x" , 

where x' and x" denote lower degree terms, we use the notation 

A(x) := ^x' i^x". 

One then sees easily that the cocycle condition (|3.ip is equivalent to 

(3.2) AB+ = {id(g)B+)A + id(g)B+{l), 

where, as in the renormalization case, with Bg'{l) = is a single vertex with assigned label S. 
The argument is then as in the original case, with the first term in the right hand side of (|3.2p 
accounting for the admissible cuts where the root vertex remains attached to the pc (T) part of an 
admissible cut of T, and the second term counts the case where the admissible cut separates the 
root vertex completely. □ 

Notice that, for the cocycle condition to hold, we need to be able to assign the same type of 
label (b, c, or r) to vertices of arbitrary valence, hence the reason for including in the set of vertex 
decorations all the fc-ary versions of composition and recursion and working with the version Hg^^ y, 
of the Hopf algebra of flow charts. 

It is possible to define other interesting grafting operators in the Hopf algebras of fiow charts, 
which, however, do not satisfy the cocycle condition. These can still be used to construct Dyson- 
Schwinger equations, but the set of solutions does not define a Hopf subalgebra. This will be 
discussed in an example in the following subsection. 

3.2. Systems of combinatorial Dyson Schwinger equations. Given a graded Hopf algebra 
H, we consider the direct product H = IlS^o^n' '^hose elements we write as infinite sums x — 
^^qX„, with Xn G Hn, endowed with the associative product and coproduct induced by those of 
H. Following |TOj, (see also [2], [30]) given a formal power series 

oo 

P{t) = ^flfci'', with ao = 1, 

fc=0 

and a Hochschild 1-cocycle i?+ on the Hopf algebra H, the associated Dyson-Schwinger equation 
is given by 

(3.3) X = B+{P{X)). 

One considers this as an equation in H and interprets the infinite sum P{X) = a^X^ as an 
element in T-L. 

One should read (j3.3p as a fixed point equation for the nonlinear transformation 

X H> B+{P{X)) 

of Ti. Solutions to the fixed point equation are elements of Ti. that exhibit a self-similarity property 
with respect to the operation of first applying X h- > P{X) (in the case of a polynomial, this 
operation heuristically replaces X by a sum of multiple copies of itself weighted by the coefficients 
of the polynomial) and then grafting them together according to the cocycle B'^. 
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In fact, the equation (13.31) has a unique solution given by an element X — Xk with Xk G 'Hk 
whose homogeneous components are determined inductively by the procedure 

n 

(3.4) Xn+l = ^ B+iXj, ■■■XjJ, 

fc=l jiH hjk=n 

starting with Xi given by i?+(l), see Lemma 2 of [5] and Proposition 2 of [lOj . 

Notice that we can still define the Dyson-Schwinger equation (|3.3I) for other types of grafting 
operators that do not necessarily satisfy the 1-cocycle condition (|3.ip . such as the examples discussed 
in the previous subsection. The existence and uniqueness of solutions, of the form p.4p is also still 
valid. Where the 1-cocycle condition is crucially used is in showing that the associative subalgebra 
of H generated by the components Xk of the solution is also a Hopf subalgebra, see Theorem 3 of 

m- 

3.2.1. Systems of Dyson-Schwinger equations. As discussed in the work of Foissy |11| . in cases like 
our "Hflo^ y/, where one has different vertex labels 6 £ {b, c, c}, one can consider more complicated 
systems of Dyson-Schwinger equations, involving the grafting operators Bf on 'Hflo^ y/. 

For 5 e {b,c, r}, we consider the grafting cocycles i?^ on the Hopf algebra W^'^^y,, as above. 
We also consider data of three non-constant formal power series in three variables X = {Xs) 

fel,fe2,fe3 

To these data we associate the system of Dyson-Schwinger equations as in [11] 

(3.5) Xs^B+{Fs{X)). 

As shown in Proposition 5 of |llj . these systems of equations also have a unique solution of 
the form Xs — "^ith the sum over all planar rooted trees with root decorated by S, with 

coefficients 



3 /•sr-^mk M 

^ _ , ,7 [1.1 = 1 PS,l)U {S) ^Pi.i... 



when 



^ -° V14 ^l,mi ^3,1 ^3, ma )■ 



3.3. Hopf subalgebras and Hopf ideals. Bergbauer and Kreimer showed that, under the as- 
sumption that i?^ is a Hochschild cocycle, solutions of Dyson-Schwinger equations determine a 
sub- Hopf algebra. Namely, they considered Dyson-Schwinger equation in a Hopf algebra % of the 
form 

oo 

(3.6) X = l + Ec„B|(A"+i). 

n=l 

One then considers the associative algebra A defined as the subalgebra of Ti. generated by the 
components Xn, with n > of the unique solution of this Dyson-Schwinger equations. In Theorem 
3 of [5], they showed that A is in fact a sub-Hopf algebra., by inductively using the cocycle condition 
p.ip to see that 

n 

(3.7) A(a;„) =En^^a;fc, where H^ = ^ ^n'-'^j.+i- 

k=0 jiH hifc + i=n-A; 
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This result was extended by Foissy to the more general form p.3p of Dyson-Schwinger equations 
in [To] and to systems of Dyson-Schwinger equations in [11]. In Theorem 4 of [10] it is shown that, 
for an equation of the form p.3|) , the subalgebra A spanned by the solutions is a Hopf subalgebra 
if and only if the formal series P{t) satisfies the differential equation 

(3.8) il-aPt)P'{t) ^aP{t), with P(0) = 1, 

for some a, /3 G K, again under the assumption that _B+ satisfies the cocycle condition (jSTTJ). 
Similarly, in 11 combinatorial conditions and conditions on the multivariable series Fi are identified 
that completely characterize when the solutions of that systems of Dyson-Schwinger equations 
generate a Hopf algebra, under the assumption that the satisfies the cocycle condition p.l|) . 

3.3.1. Ideals of Dyson-Schwinger solutions. It is natural to consider not only the associative sub- 
algebra A spanned by the components of the solutions of Dyson-Schwinger equations, but also the 
ideal X spanned by the (xn)n>\- As in the case of the subalgebra A^ it is also natural to ask for 
conditions that will ensure that it is a Hopf subalgebra, so with the ideal X of solutions it is natural 
to ask for conditions ensuring that it is a Hopf ideal, so that it makes sense to define a quotient 
Hopf algebra, playing the role of the ring of functions on the variety of solutions. 

As an example, we consider here the Dyson-Schwinger equation p.3p . with i?+ = the grafting 
operator to a vertex labelled with the bracketing operation. Assuming that parallel computations 
that do not feed into each other can be performed in arbitrary order, we see that this grafting oper- 
ator descends to the commutative quotient 7^flow,v' of the noncommutative Hopf algebra 'W^^ y,. 
Let X — X)n>i be the unique solution of p.3|) . obtained as in p.4|) . with x\ = B'^{X) = «[, a 
single vertex labeled with the b operation. We define X to be the ideal in the commutative algebra 
'Hflow,v' generated by the a;„ with n > 1. 

The same conditions that ensures that is a Hopf subalgebra also ensure that Z is a Hopf ideal. 
We verify it in the case of (j3.6p . The more general case of [10 is similar. 

Lemma 3.2. The ideal X generated by the components Xn with n > 1 of the solution of p.6p is a 
Hopf ideal. 

Proof. Elements of the ideal X in H — Hfiowy are finite sums X]m=i iT-mXk^, with h,n E % and 
Xk the coordinates of the unique solution of p.6p . The condition that I is a Hopf ideal is that 
A (J) C I (gi H © H (g) I. Using the formula (13.71) for the coproduct of the elements ccfe, we see that 
we obtain a sum of terms of which the primitive part l®Xk+Xk®l \s, vn.T-L®X(BX®T-L and all the 
other terms K{xk) are uyX®X. Then the coproducts A{hmXk^) will also be in T-L^X^X^I-L. □ 

Notice that we can also work with the noncommutative Hopf algebra 'H^'^^ y, and the two sided 
ideal with elements of the form hmXk^tm, with hm,(-m G ^how V' ^'^^ considerations 
that follow (see ij3. 41 below) it is more natural to work with the commutative quotient. 

3.4. Yanofsky's Galois theory of algorithms. As we have seen above, one can pass to commu- 
tative quotients Haow of the noncommutative Hopf algebras of flow charts. The meaning of passing 
to the commutative quotient Haow is best expressed as follows: one can think of monomials in "HJo^ 
as diagrams for parallel computations of a certain number of outputs (one for each planar tree in 
the forest). Imposing commutativity then corresponds to the reasonable expectation that if one 
has two parallel computations that do not feed one into the other, then computing them in either 
order will not affect the result. However, when one introduces operations that graft the output of 
one three to the input of another, with c and r labels at the vertices, imposing the commutativity 
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relation is no longer appropriate, as the result of grafting would not be well defined as a planar 
tree. In the case of grafting with c label, we can still pass to the commutative quotient, as argued 
above. 

3.4.1. Commutative Hopf algebras and affine group schemes. An advantage of working with com- 
mutative Hopf algebras is that they are dual to afhne group schemes. Thus, the Hopf algebra Hhow 
determines an affine group scheme Gnow such that, for any commutative K-algebra A, one obtains 
a group G'flow(^) = Hom(Hfiow, the set of homomorphisms of K-algebras, with the product 

-k(j)2){x) = (X) 02, A(a;)) dual to the coproduct and the inverse determined by the antipode. 

3.4.2. Galois groups of algorithms. The idea elaborated by Yanofsky in the recent work |29| . of a 
Galois theory of algorithms, aims at considering all possible "reasonably sets of relations" on the 
Hopf algebra Ti^^^ of flow charts that would correspond to "implementing the same algorithm" , 
viewing "algorithms" as an intermediate level between the labelled planar rooted trees and the 
primitive recursive functions they compute. 

Reasonable relations between the planar rooted trees of TL^^^ correspond therefore to passing to 
quotient Hopf algebras. Each such quotient is determined by a Hopf ideal I in TL^^^, the kernel of 
the map to the quotient Hopf algebra Hi. 

In the commutative setting T^aow one can view, dually, the quotient Hopf algebras Hi as sub- 
group schemes G/ of the affine group scheme of characters Gflow The group schemes Gi obtained 
in this way are the basic objects of interest in Yanofsky's Galois theory of algorithms. 

3.4.3. Hopf ideals and Galois groups. In light of Yanofsky's point of view, considering ideals in Haow 
generated by solutions of Dyson-Schwinger equations, and the resulting quotient Hopf algebra, 
rather than sub-Hopf algebras, corresponds to considering the "ring of function" of the variety of 
solutions. Namely the associated group scheme G/ consists of those symmetries in Gfiow that are 
compatible with the self-similarity structure described by the Dyson-Schwinger equation. 

Proposition 3.3. A Dyson-Schwinger equation p.6|) in the commutative Hopf algebra Hflow de- 
termines a Galois group Gp C Gflow 7 in the sense of Yanofsky 29 . 

The result follows directly from Lemma [3.2l A similar result holds for Dyson-Schwinger equations 

(13:3)) with dsii). 

3.5. A non-Hopf example with binary trees. In the setting of Dyson-Schwinger equations 
described above, it is crucial that we consider the planar rooted trees in Hf[^^ y, without external 
edges. In fact, this allows us to define the grafting operators B'^ as grafting any number of trees to 
the same root consisting of a single vertex decorated by S. If we reintroduce external edges, then 
the graphs consisting of a single vertex but with different numbers of incoming edges are viewed as 
different generators in the Hopf algebra and they no longer define a cocycle i?^. 

To see this more explicitly, we focus on the case of binary trees, keeping track of external edges. 
Consider the following primitive elements in the Hopf algebra 'Hg^^y: 



(3.9) 
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For each of these graphs r one can define a grafting operation i?+ by setting the result equal 
to zero on monomials T that cannot be grafted to r. This _B+ clearly does not satisfy the cocycle 
condition because one would have to have B^{1) = 0, but then the expressions Ai3+ and A, 
when applied to elements T = n 11x2 that are not in the kernel, would differ by a term corresponding 
to the cut that separates r from T. On the other hand, if B^{1) = r as this case would require, 
then the cocycle condition fails on elements in the kernel of 5+ . 

Another way of defining a grafting operator with r one of the trees (|3.9p is to define, for a 
monomial T = ti 11 • • ■ 11 r„ the grafting B^[T) to be obtained by just gluing one of the outputs 
of T (say, the first one) to one of the two inputs of r (say, the first one). We use this choice in 
the explicit example below. Again, however, one sees that the grafting defined in this way does 
not satisfy the cocycle condition. In fact, consider the case of a monomial T consisting of a forest 
with more than one component. Then only one of the components of T is glued to r by B^ . This 
means that, in the computation of the term AB+ there is not just one but several admissible cuts 
that separate r from T, namely all the admissible cuts that include a cut of the edge used for 
grafting together with an arbitrary cut of the remaining components not connected to t. Thus, the 
difference between Ai?+ and [id ® B^)K does not consist only of the term id ® Bf[l) — id®T. 

However, one can still consider a Dyson-Schwinger equation of the form 

X = l + Y,B+{X^), 

s 

where the sum is on the three labels (5 G {b, c, r} of the vertex of the tree ts as in p.9p . This admits 
an explicit solution, as in the general case above, which are of the form 

n 

Xn+1 = X! X! {^kXn-k)- 
k=0 S 

In Hfiow V' one can see that the first few terms are of the form xq = I and 



Xi = 
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+ 



+2 



+ 



+2 



One can continue with the successive coefficients, such as = B^{x\+2x2Xq) , using the facts 
that for a vertex iabeled by recursion, the order of inputs always matters, as well as for composition, 
while we can assume that order of inputs does not matter for bracketing as a consequence of the 
nature of parallel computation. 

This example shows that, if we want to work with binary trees with external edges, which is 
a very natural choice from the point of view of computation theory, and we want to obtain a 
grafting operation that is non-trivial on monomials of arbitrary degree, we would have to extend 
the grafting to include, besides the primitive two-valent trees (|3.9|) . also a choice of binary trees with 
larger numbers of input flags. This cannot be accommodated in the usual setting of combinatorial 
Dyson-Schwinger equations based on Hochschild cocycles in Hopf algebras. However, we shall see 
below that such generalizations exist naturally when one reformulates Dyson-Schwinger equations 
in the operadic setting. 



4. Operadic viewpoint 

Instead of using a Hopf algebraic setup to describe flow charts and primitive recursive functions, 
one can also adopt an operadic viewpoint, as suggested by Manin in |24| . 

Let us consider, as above, planar rooted trees r with vertex labeling by the elementary operations 
of type {b, c,r}, oriented so that each tree has a single outgoing flag attached to the root vertex 
and a certain number of incoming flags. 

Definition 4.1. We define the operad of flow charts Oflow.v setting 0{n) to the the 'K-vector 
space spanned by labelled planar rooted trees with n incoming flags and the operad composition 
operations 

OQ : 0{n) ® 0{mi) 0{m,i) 0{mi H h m„) 

are given on generators r ti (g) • • • iX) t„ by grafting the output flag of the tree ti to the i-th input 
flag of the tree r. 
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As pointed out in 2 , Dyson-Schwinger equations admit a natural operadic interpretation. 
Namely, given a formal series P{t) = 1 + J^'kLi '^^t'' and a collection /3 = (/3„) with /?„ G 0{n), we 
consider the equation 

(4.1) X = (3{P{X)), 

where X = J^k is a formal sum with Xk G 0{k). In the right hand side we have P{P{X))i = 
1 + (3i o xi, where 1 is the identity in 0(1), and for > 2 

n 

(4.2) /3(P(X))„=^ ^ ak (3ko{xj,(^---®Xj,), 

k=l jiH hjk=n 

with Xj^ ® • • • ® Xj^ e 0(ji) (8) • • • (8) 0{jk), so that the composition (xj^ lE) ■ ■ ■ (S) Xj^) e 0(71), 

since ji H h jfe = n. 

Proposition 4.2. For O ~ Ofiow.v the operad of flow charts, if ai/3i ^ 1 G the operadic 

Dyson-Schwinger equation (14. ip /las a unique solution X G nn>i ^l'^) g^'ven inductively by 

n+l 

(4.3) (1 - ai/3i) o a;„+i = ^ ^ Ofe ^fe o (x^^ 8) • • • 8) Xj J. 

fe=2 j"i+---j"fc=n+l 

Proof. First observe that in the operad of flow charts 0(1) is one-dimensional, since it is spanned by 
the trees with a single vertex, one incoming and one outgoing flag, with the operations b, c or r at the 
vertex, but in the case of a single input, all of these operations are the identity, so that elements of 
0(1) are scalar multiples of the identity map, viewed as operations on primitive recursive functions. 
Thus, the element (3i G 0(1) is a scalar multiple of the identity and, if we assume that ai/?i ^ 1 in 
0(1), then 1— ai/3i is invertible. Since we have /3{P{X))i = l+/3ioxi as the term in 0(1) in the right 
hand side of (|4.ip . we obtain f 1— ai 0i )xi = 1, which fixes xi = {l—ai(3i)~^ S 0(1). At the next step 
we have, from (|4.2p . X2 — ai/3ioa;2 + a2/32°(a;i8)Xi), which gives X2 = (1 — ai/3i)^"'^(a2/32°(2^i'82;i)), 
with xi as above. At the (n-l- l)-st step, Xn+i is then determined uniquely by (j4.3p in terms of the 
coefficients ak and the elements /3k for fc=l,...,n-|-l. □ 

Let then Op^pin) denote the K-linear span of all the compositions Xk o {xj-^ 8) ■ ■ • 81 Xj^^), for 
k = 1, . . . , n and ji + ■ ■ ■ + jk = n, with Xk the coordinates of the solution X = Xk of the 
Dyson-Schwinger equation determined by /3 = (/3„) and P{t). The Op^p{n) form a sub-operad 
under the composition maps induced from 0(n). 

For instance, for a choice of one of the operations 5 € {b,c,t}, we can take oi 7^ 1 and the 
element Pk given by the tree with a single vertex marked by (5, with k incoming flags and one 
outgoing flag. This choice gives an operadic reformulation of the Hopf-theoretic Dyson-Schwinger 
equations with the cocycles B'^ . However, we can now consider also more general operadic Dyson- 
Schwinger equations for different choices oi (3 — (/3„) that do not correspond to Hochschild cocycles 
in the Hopf algebra setting. 

Restricting to the sub-operad O^^p(n) generated by solutions of the Dyson-Schwinger equation 
should be regarded as considering those operations (formal combinations of flow charts) that satisfy 
a self-similarity property with respect to the transformation X (3{P{X)). 
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4.0.1. The case of binary trees. In particular, we can now revisit the example of the binary trees in 
this operadic setting. Let C^o^ v(^) K-linear space spanned by the planar binary trees with 

vertex labehng by the binary operations b — b2, c. These form a sub-operad of Oflow.v with the 
induced composition maps. 

For a fixed S £ {b, c, c} let f3i^s be the binary tree with a single vertex, two incoming and one 
outgoing flag, as in p.9|) . We then consider a choice /9„ of a binary tree with n input flags, for 
each n > 1, with a fixed choice of vertex labelings in {b,x, c}, not necessarily all of the same type. 
Each such choice oi f3 — {I3n), together with a choice of P{t), now determines an operadic Dyson- 
Schwinger equation of the form (j4.1l) . in the binary setting, with a corresponding solution as in 
Proposition 14.21 

4.0.2. Systems of Dyson-Schwinger equations in operads. By analogy to what happens in Hopf 
algebras, 11 , one can pass from the case of a single Dyson-Schwinger equation to systems of 
Dyson-Schwinger equations. In the operadic setting a system of Dyson-Schwinger equations would 
be determined by the data of formal series Fs{ti,t2, is) = J2ki fea ^^ki k2 fe3^i^*2^*3^ ^^"^ ^f elements 
= (/Si*^^), where for 5 G {b, c, r}, the elements Pn'' G 0{n) are realized by rooted trees with root 
vertex marked by 5. The operadic systems of Dyson-Schwinger equations are then of the form 

Xs=P^'\Fs{X)). 

The explicit form of the solutions is more cumbersome than in the case of a single equation in 
Proposition 14.21 above, and it follows a pattern similar to the derivation of the explicit solutions of 
Hopf theoretic systems of Dyson-Schwinger equations p.5|) as in [11] . These systems detect more 
elaborate forms of self-similarity in the operad of flow charts, involving different families of grafting 
operations simultaneously. 

4.0.3. Operads and Properads. As mentioned in [24], one can extend the Hopf algebra of flow charts 
by considering, instead of labelled embedded rooted trees, more general labelled embedded acyclic 
graphs, namely graphs endowed with an acyclic orientation. Correspondingly, the operad of flow 
charts would be replaced by a properad, where the compositions extend from grafting output and 
input flags of trees to grafting outputs and inputs of acyclic graphs. 

The notion of properad was introduced in [27) as an intermediate notion between operad and 
prop. Namely, a properad parameterizes operations with varying numbers of inputs and outputs 
that are labelled by connected acyclic graphs, while the operad case uses trees (varying number of 
inputs and a single output) and props allow for disconnected graphs. 

More precisely, we consider the properad T^sow.v where Virn, n) is the K-vector space spanned 
by planar connected directed (acyclic) graphs with m incoming flags and n outgoing flags, with 
vertices decorated by operation that include the elementary b, c, r, with m inputs and one output, 
but also include now a chosen set of additional operations with m inputs and n outputs. As we 
have seen in the case of extending the operations c and c from binary to m-ary inputs, in this case 
also such additional operations with m inputs and n outputs acting on primitive recursive functions 
can be decomposed in terms of the elementary operations: the resulting operations associated to 
vertices with m incoming and n outoing flags can be regarded as "macros" , in the sense discussed 
in §2 of [29j . The properad has composition operations 



r{m,n)®V{ji,ki) ®---®V{jt,ke) V{ji + ---+je,n), for ki + ■ ■ ■ + h = m. 
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4.0.4. Dyson-Schwinger equations in properads. We can define Dyson-Schwinger equations for pr- 
operads, extending tlie operadic case considered above. Namely we consider again a clioice of a 
formal power series P{t) = l+J2k ^kt'^ ^^'^ ^ collection P — {(im,n) of elements with Pm.n & Vim, n). 
We then consider the equation 

(4.4) X - P{P{X)) 

as in (j4.ip , where now the right hand side has (m, n)-componcnt in V{rn^ n) given by 

in 

(4.5) /3(P(X))„,„-^afc 

k=l - jfc=m 

In order to construct solutions of the properad Dyson-Schwinger equations, we introduce trans- 
formations A„ = A„(a, ji) with 

A„(a,/3) : ©fc^iP(n,fc) -> ©^^iP(n,fc), with A„(a, /3),j = a^^j- 



Theorem 4.3. If for all n > 1 the transformation / — A„(a, /3) is invertihle, with I the identity on 
©fc=i^("i fc), then the properad Dyson-Schwinger equation (|4.4p admits a unique solution, given by 
Xi^i = Aj~^, and for m < n by 

m 

(4.6) Xm.n = y^^a-kl^k, 



( k 



k=l 



E E 

= 1 3lH hjf=™ 

iiH hi, = k 



(8) ■ • • (8) 



while the remaining components Xm.n with m > n are determined by 
(4.7) r„(:r) = (/-A„)-iA„y(")(x), 

where y„(x)* = (x„,i, . . .,Xn,n) and F(")(a;)* = (l^(")(a;)j)j=i,...,„ 



^^"H-).=E E 



Xr,.,. 



k=2 ri + .-.+rfc=„ 
= l+--- + = fc=3 



Proof. By (j4.5[) we see that we have xi^i = 1 + ai/3i.i o xi^i, which can be solved for xi^i provided 
the invertibility of 1 — Ai = 1 — ai/3i_i holds. The components Xm,n with m < n are determined 
by (j4.6p . where the right hand side only involves components Xr,s with r < m and s < n. The 
equations for the remaining entries Xm,n with n < m break into systems of the form 



/ X' 



{I - A™) 



m,l 



= E(A")dE E 



■ (g) Xr 



fc = 2 '■! + .. . + rj^=™ 
= l + --- + = fc=J 



Provided that the invertibility condition for the transformation I — Am holds, these systems provide 
unique values for all the components Xm,n with n < m, as a function of the components Xr^s with 
r < m, that we have already been determined. □ 

Notice that one does not expect the invertibility of / — A„(a, /3) to hold in general, as A„(a,/3) 
involves compositions with the elements Pj^i. The simplest example in which one has invertibility 
can be obtained by taking l3m,n = for m ^ n and / — an invertible transformation in V(n, n). 
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As in the operad case, we can define Vp^p{m,n) to be the K-vector spaced spanned by all the 
properad compositions 

with ji + ■ ■ ■ + ji = m and ii + ■ ■ ■ ii = k and where all the Xm,n are components of the solution to 
the properad Dyson-Schwinger equation determined by the choice of P{t) and of /3 = (/3m, n), with 
the invertibility condition of Theorem 14.31 The Vi3^p{m,n) form a sub-properad of "Pflow.v with 
the induced composition. 



5. RENORMALIZATION of the halting PROBLEM 

In [22], Manin adapted the formalism of perturbative renormalization to the halting problem in 
the theory of computation. To summarize Manin's approach, the basic idea is that one considers 
non-computable functions as an analog of divergent Feynman integrals. Renormalization consists 
of a procedure of extraction of finite values from divergences and in this setting it corresponds to 
extracting a "computable part" from a non-computable function. 

5.1. The halting problem: regularization and renormalization. We describe here a setting 
for renormalization of the halting problem which is minor elaboration on the procedure described 
in §3 of dg. 

5.1.1. Partial recursive functions and Hopf algebras. We first pass from the class of primitive re- 
cursive functions considered above to the larger class of partial recursive functions. These have a 
similar presentation as the one recalled in 321 above, with the same basic functions (successor, con- 
stants, and projections) and, in addition to the three elementary operations c, b,r of composition, 
bracketing and recursion, an additional ^ operation that assigns to an input function / : N"^^ — N 
assigns an output 

/i : N" N, h(xi,. ..,Xn)= min{a;„+i | /(xi, . . . ,x„+i) = 1}, 

with domain 

2?(/i) = {{xi,...,Xn) \ 3xn+i > 1 : f{xi, . . . ,Xn+i) = 1, with (xi, . . . , a;„, fc) e 'D{f),\/k < x„+i}. 

According to Church's thesis the partial recursive functions obtained in this way are exactly the 
semi-computable functions, namely those for which there exists a program that, for all x £ T^{f), 
computes f{x) and which either computes zero or runs for an infinite time when x ^ ^5(/), see §V 
of [2n| • Notice that it is the presence of the additional operation /i that allows for the construction, 
starting from the basic functions, of functions that are only partially defined and semi-computable. 

Correspondingly, we enrich the Hopf algebras of flow charts by additional vertex decorations 
by the /i operations, extended to vertices of valence greater than two by a combination with the 
bracketing operation. We shall still use the notations 'H^'^^ ^, and H^^^ y, , for the resulting Hopf 
algebras with admissible decorations by partial recursive functions or with vertex decorations, 
respectively, and 'Hfiow,K' and Hflow.v for their commutative quotients. 
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5.1.2. Algebraic renormalization. The basic algebraic formalism for perturbative renormalization, 
in its formulation as given in [9] (see also §5 of [25]), consists of an algebraic Feynman rule 

(5.1) (jj-.n^B 

which is a morphism of commutative algebras from a commutative Hopf algebra Ti to a Rota-Baxter 
algebra B of weight — 1 . 

Recall that a commutative Rota-Baxter algebra of weight A is a commutative, associative algebra 
B endowed with a linear operator T on it that satisfies the identity 

(5.2) T(x)T(y) = T{xT{y)) + T{T{x)v) + A T{xv). 

The field of convergent Laurent series in one variable with the projection T onto the polar part 
is an example of Rota-Baxter algebra of weight A = — 1, which is widely used in renormalization 
theory, 'S], [5]. The Rota-Baxter operator T of weight —1 determines a splitting of B into two 
commutative unital subalgebras, S+ = (1 — T)B and B- given by TB with a unit adjoined. The 
fact that B± are indeed algebras and not just vector spaces follows from the Rota-Baxter identity 
(1521). 

The BPHZ renormalization of the algebraic Feynman rule (jS.ip is then obtained by showing the 
existence of a multiplicative factorization 

(/)=((?!)_ o 5) with (t>±:n^B±, 

obtained inductively via the BPHZ preparation formula (see [3], [9]) 

(5.3) 0_(x) = -r(0(x)+^0_(x')0(x")), and 0+(x) - (1 - r)(0(x) + ^ 0_ (a;')0(x")), 

where a; e H, with A(a;) — l(x)x + x(x)l + '^x'<Si x" . The factorization is unique if normalized 
by €- o — e, where e_ : B- — > C is the augmentation map and e is the counit of H. Again, 
it is the Rota-Baxter property of {B, T) that ensures that 0± obtained in this way are algebra 
homomorphisms. The homomorphism : % B+ is the renormalized algebraic Feynman rule, 
while (ji-T-i ^ B- is the counterterm, the divergence. 

5.1.3. An algebraic Feynman rule for the halting problem. We let B be the algebra of functions 
<j) . pjm _^ M[D) from N™, for some k, to the algebra J^{D) of analytic functions in the unit 
disk D = {z e C : |z|<l}. Since the functions $ can have different domains N*^, the sums and 
products are defined by extending them to functions ^{ki, k2, . . . , fcn, • . .) from a common domain 
N°°, which depend on only finitely many variables, $ = $ o tt^, with tt^ : N°° — >■ N™. The Rota- 
Baxter operator T on S is given by componentwise (for each k) projection onto the polar part at 
z = 1. Thus, B+ consists of functions from N°° (depending only on finitely many ki) to Ai{D)^, 
the algebra of meromorphic functions in D that extend continuously at z = 1 , while B- consists of 
functions that are either constant or have a pole at z = 1. 

We define algebraic Feynman rules for the Hopf algebras Haow.TC' and Haow.v- In the case of 
"Hflow.K', where both vertices and flags in the trees are labeled by admissible labelings, there is a 
well defined partial recursive function / associated to each labelled rooted tree r in 'Hiiow,n' , namely 
the function / labeling the unique outgoing flag of r. Let us first assume that the partial recursive 
function / has range in N. We extend this / : N"* — ^ N to a function / : N™ — ^ Z>o that computes 
f{x) at X G !'(/) and takes value at x ^ 2?(/). 
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We assign to a tree r that computes / an element in ^rik, z) in defined (as in §3 of [22]) by 
(5.4) $.(fc,z)^$(fc,/,z):^^ 

The resulting function ^rihr) G A^(-D) has a pole at z = 1 iff fc ^ T^if)- We extend this 
definition mutiplicatively to cases where the output of the flow chart r is a partial recursive function 
/ : N*^ ^ N^, by setting $r(fc; z) — 0^=1 ^{h^ fji^)- In tbis way, the resulting function has a pole at 
z = 1 if at least one of the fj does, that is, iff fc ^ Z^(/i) H • • ■r\V{fe) — !?(/), consistently with the 
above. We then extend the definition multiplicatively to the case of a monomial ti • • • t„ in Hflow.-R.' 
and additively to linear combinations of monomials. Thus, we obtain an algebra homomorphism 

We write $(A:, /, z) :— $/(fc)(z) for the function ^f{k) G M{D) oi z € D, associated to a given 
k € N°° and depending on the output / of a given flow chart, with / a partial recursive functions. 
If we denote by r the decorated planar rooted tree describing the flowchart, we equivalently write 

C&(fc,T, z) or $fe,r(2)- 

The renormalization of the halting problem described by Manin in f2? then consists of applying 
the BPHZ factorization ()5.3p to these algebraic Feynman rules to consistently extract a computable 
part from the halting problem for partial recursive functions that accounts for "subdivergences" 
created by subroutines within the flowcharts. 

When we apply the BPHZ procedure to <i> defined as in (|5.4|) . we obtain 



(5.5) <f _(fc, fr, Z) = -rmk, fr, Z) + 2^ <i>-{k, U^^r),zMk, fpc{r),z)). 

C 

Because f = fr is the output function computed by the flow chart, that is, the label of the outgoing 
flag of r, we have fpdr) = fr, since these trees have the same root and outgoing flag with the same 
label. Thus, the expression above is simply of the form 



(5.6) 



^-{L fr, Z) = -T (^{k, fr, Z)(l + $_(fc, ^ f^dr),z))^ 



This has the effect of considering not only the pole at z = f that is caused by fr itself, but also 
those that come from subdivergences created by all the partial recursive functions f-jrc^r) computed 
in the intermediate steps of the computation performed by the flow chart t, which are the outputs 
of the "pruned parts" ■nci'r) of the admissible cuts on r. 

5.1.4. The case of T-Lfio^^\;> . A variation on the construction above would be to define an algebraic 
Feynman rule on 'Hflow.v'i where the flags are not labelled and only vertices are. In this case, 
in order to define something like (15.41) . we need a way of assigning inputs to the flow charts. A 
reasonable choice is to use the basic functions as input. We then define the algebraic Feynman rule 

as 

(5.7) Hk,r,z):=l[mJr,a,z), 

a 

where the product is over a ranging over the set of all possible functions tr : EI"^(t) — ?• {s,c, tt} 
that label the incoming external edges of r by basic functions (successor, constant, or projection). 
The partial recursive function fr^a is the output of the flow chart given by the tree r with inputs 
assigned by a. Each $(fc, fr,a, z) is computed as in (|5.4|) . and ()5.7|) is extended to arbitrary elements 
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of the Hopf algebra as before. In this setting, the function ^rik, z) has a pole at z = 1 iff there is 
at least a choice of a basic input a for which k ^ ^{fr.cr)- 

With this setting, the BPHZ formula for the algebraic Feynman rule (|5.7p becomes more inter- 
esting than in the case of (|5.6p . We have again (|5.5I) . in the form 

(5.8) (fc, r, z) - -T($(fc, r, z) + ^ (fc, ^c(r), z)<i>(fc, Pc(r), z)), 

c 

with $(A:, pc'('''), z) and $(A;, 7rc(T), z) again computed as in (|5.7p . but this time it is no longer 
true that fr = fpc(T) because on the tree pc{t) we are using the new inputs given by basic 
functions and not the input coming from the output of 7rc(T). The divergence (jS.BI) here combines 
the divergences of the fr^a with divergences coming from partial recursive functions /7rc(r),a and 
fp^.(r),<j that flowcharts ttc{t) and pc{t) compute starting from inputs of basic functions. 

5.2. Dyson Schwinger equations for the halting problem. If we work with the Feynman 
rule described above on the Hopf algebra Hhow,!/" with only vertex decorations, then we can read- 
ily consider Dyson-Schwinger equations as discussed in iJ3] above. However, if we work with the 
definition of the algebraic Feynman rule $ described in (|5.4[) on the Hopf algebra 'Hfiow,TZ' of flow 
charts with flag decorations, we need to extend the grafting operators appropriately that are needed 
to define Dyson-Schwinger equations from Haowy to Hflow.v- This can be done by assigning as 
output the empty function / = 0, mapped to the constant function ^{k, 0, z) = 1, as the output of 
any tree r that is obtained through a grafting where the flag labels don't match. 

We can then consider a Dyson-Schwinger equation of the form p. 31) . with i?+ = X^i-^/' with 
S € {b,c,r, /i} or a more general system of Dyson-Schwinger equations p. 51) . with the condition 
that the components Xm of the unique solution span a Hopf subalgebra of T^aow.TC'- In such case, 
we can restrict the algebraic Feynman rule $ constructed above to this Hopf subalgebra and still 
perform the BPHZ renormalization. The divergences $-(&, fx^i , z) (or '^-{k, fx,„,a-, z) in the case 
of Hflow.v) measure the "amount of non-computability" that can be produced by flow charts that 
satisfy the self-similarity property with respect to the operations {b, c, r, /i} specified by the Dyson- 
Schwinger equation. 
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